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Various types of active galactic nuclei (AGN) are briefly discussed, with an emphasis on the 
theory of recent X-ray and 7 -ray observations of the subclass, Seyfert 1 galaxies. The large 
radiation power from AGN is thought to originate from gravitational power released by matter 
accreting onto a supermassive black hole. The physical mechanisms responsible for the energy 
release and the geometry of the gaseous components are still uncertain in spite of three decades 
of observational and theoretical studies. The recent X-ray and 7 -ray observations, however, start 
to provide useful constraints on the models. This kind of interpretation of the observations is 
possible due to theoretical developments during the last decade of radiative transfer of X-rays in 
both optically thick and thin media of various geometries. Particular attention is paid to various 
accretion disk-corona models. Recent work on exact radiative transfer in such geometries are 
reviewed. 


1. Theoretical history 

The underlying framework for almost all efforts to understand active galactic nuclei 
(AGN) is the accretion disk picture described in the classical papers by Novikov & 
Thorne (1973) and Shakura fc Sunyaev (1973)| . This original picture was partly inspired 
by the extreme optical AGN luminosities. Here, effectively optically thick, rather cold 
matter forms a geometrically thin, differentially rotating Keplerian disk around a super- 
massive black hole. The differential motion causes viscous dissipation of gravitational 
binding energy resulting in outward transportation of angular momentum and inward 
transportation of matter. The dissipated energy diffuses vertically and emerges as black 
body radiation, mostly in the optical-UV spectral range for the case of AGN. 

Observations have also been the driving force in the discovery of two other solution 
branches. 

a) Hard X-rays from the galactic black hole candidate, Gyg X-1, as well as the discovery 
of strong X-ray emission from most A GN led to the need for a hot accretion disk solution. 
^hapiro, Lightman fc Eardley (1976) (SLE) found a hot, effectively optically thin, rather 
geometrically thin solution branch, where the ions and the electrons are in energy balance, 
with the ions being heated by dissipation and cooled through Goulomb exchange, leading 
to ion temperatures of order 10^^ - 10^^ K. The efficient cooling of electrons through a 
variety of mechanisms above 10® K, leads to electron temperatures being locked around 
10® - 10^® K. The SLE-solution is thermally unstable. 

b) Observations of Gyg X-1 and of radio galaxie s led to two independent discover¬ 
ies of the third solution branch. Ichimaru (1977)| developed a model to explain the 
soft high and hard low state of Cyg X-1. The soft high state is due to the disk be¬ 
ing in the optically thick cold state, and the hard low state occurs when the disk 
develops into a very hot, optically thin state. This solution branch is similar to the 
SLE-solution, except that now the ions are not in local energy balance. It was found 
that if the ions were sufficiently hot they would not cool on an inflow time scale, but 
rather the ions would heat up both by adiabatic compression and viscous dissipation 
and would carry most of that energy with them into the black hole. Only a small frac¬ 
tion would be transferred to the electrons, so the efficiency of the accretion is much less 
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than the normal ~ 10%. As the ion temperature was found to be close to virial, these 
disks are geometrically thick. Just as for the SLE-case above, the electrons decouple 
to be locked at 10® - 10^® K. The flows resemble the quasi-spherical dissipative flows 


studied by, e.g., Meszaroz (1975) 

and Maraschi, Roasio & Treves (1982), Independently, 

Rees et al. (1982) and 

Phinney (1983) proposed that a similar inefficient accretion disk 


solution is responsible for the low nuclear luminosities in radio galaxies with large radio 
lobes and thus quite massive black holes. These geometrically thick disks were named ion 
tori. Rees et al. (1982)| specified more clearly than Ichimaru (1977)| the critical accretion 


rate above which the flow is dense enough for the ions to cool on an inflow time scale 
and the ion tori-branch would not exist. Further considerations of ion tori were made 


by Begelman, Sikora fc Rees (1987) . Ichimaru (1977) emphasized that his version of ion 
tori is thermally stable. 

The ion-tori branch has been extensively studied and applied over the last two years 
(1995-1996) with more than 30 papers by Narayan and co-workers, Abramowicz and co¬ 
workers, as well as many others. These studies confirm and extend the original results, 
although some papers do not quote or recognize the original results. New terminology has 
been introduced based on an Eularian viewpoint rather than a Lagrangian. Instead of the 
ions not cooling, it is said that the local volume is ’’advectively cooled” due to the ions 
carrying away their energy. The ion tori are therefore renamed as advaction-dominated 
disks. 


The three accretion disk branches were ’’unified” by Chen et al. (1996) and Bjornsson 
et al. (1996), who determined where in parameter space the branches merges. Parallel to 
these developments there have been several other lines of research regarding the physics 
needed for realistic modelling of accretion flows. These areas include radiative processes 
in hot plasmas, radiation transfer in hot plasmas, radiation transfer of X-rays in cold 
plasma, MHD in differentially rotating gas, the origin of disk viscosity, magnetic flares, 
gas or MHD-simulations of flows, and so on. Some of these research lines have already 
merged. Here we focus on the areas most important for interpreting the high energy 
spectra of AGN. 

While the classical rates for bremsstrahlung, cyclo/synchrotron radiation, and Comp¬ 
ton scattering in the nonrelativistic and relativistic limits were sufficient in the early disk 
models, the electron temperatures of 10® - 10^® K indicated by both observations and 
theory required the calculation of transrelativistic rates of the above processes as well 
as for pair processes that becomes important at these temperatures. Rate calculations 
as well as exploring the properties of pair and energy balance in hot plasma clouds were 
done in the 1980s by Lightman, Svensson, Zdziarski and others. Methods to incorporate 
these sometimes complex hot plasma solutions into the simplest accretion disk models 
were developed by Bjornsson & Svensson (1991, 1992). The Compton s cattering kernel 


probably received its definite treatment in Nagirner fc Poutanen (1994) 


To obtain approximate spectra, the Kompaneets equation with relativistic correc¬ 
tions and with a simple escape probability replacing the radiative transfer is sufficient 
(Lightman & Zdziarski 1987). However, if one want to obtain constraints on the ge¬ 
ometry from detailed observed spectra, then methods to obtain exact radiative trans- 
fer/Comptonization solutions in different accretion geometries must be developed. Here 
again, Igor Novikov had great influence as he and his group stimulated a particle physicist 
and expert on Monte Carlo simulations, Boris Stern, to develop a Monte Carlo code able 
to treat all the high energy radiative processes occurring in AGN. The first results were 
reported in Kardashev, Novikov fc Stern (198'^ and Novikov & Stern (1986), and an ad¬ 
vanced version of the code was finally documented in Stern et al. (1995a). Faster meth¬ 
ods were developed by Haardt (1993)| (approximate treatment of the Compton scattering) 
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and Poutanen & Svensson (1996) (exact treatment) who solved the radiative transfer for 
each scattering order separately (the iterative scattering method). These codes are fast 
enough to be implemented in XSPEC, the standard X-ray spectral fitting package, and 
exactly computed model spectra can now be used when interpreting the observations. 
So far these methods have been used to study radiative transfer in two-phase media con¬ 
sisting of cold and hot gas with simple geometries, but they have yet to be integrated in 
the standard accretion disk models. 

UV and X-ray observations indicate the coexistence if both cold and hot matter in 
AGN. Then one must be able to compute the reprocessing (both absorption, reflection, 
and transmission) of X-rays by the cold matter. The first approximate considerations of 


Compton reflection by |Guilbert fc Rees (1988)| and [White, Lightman, fc Zdziarski (1988) 


have now led to quite accurate methods developed by [Magdziarz k, Zdziarski (1995)| and 


Poutanen, Nagendra fc Svensson (1996) Such Compton reflection forms an important 


ingredient in the radiative transfer calculations in two-phase media. 


2. Observational history 

The exploration of the spectra of AGN started some 30 years ago. It is mainly due 
to the developments in instrumental and satellite technology that we now in the 1990s 
start having a full broad band picture of these spectra. In 1960s groundbased observa¬ 
tions started exploring the optical and radio properties. The advent of X-ray astronomy 
allowed the 2-10 keV properties to be explored in the 1970s. The theoretical modelling 
then consisted of connecting the radio, optical, and the X-ray data points with one or 
two power laws and discussing the resulting fits within nonthermal models where rela¬ 
tivistic power law electrons generate both power law synchrotron radiation and power 
law Compton-scattered radiation. The infrared spectral gap was covered in the 1980s. 
First after the launch in April 1991 of the Compton Gamma Ray Observatory (CGRO) 
covering the gamma ray spectral range did we get a full broad band view of AGN spectra 
and thus reliable estimates of their total lumino sities and the spectral rang es where most 
of the luminosity was emitted. As described in Dermer fc Gehrels (1995) , AGN fell into 
two distinct classes, the y-loud and the y-weak AGN. The y-loud AGN are the lumi¬ 
nous blazars where we are looking down into a relativistic jet emerging from the AGN. 
Here, the power emerges in two broad band humps, one due to nonthermal synchrotron 
emission peaking in the IR-optical spectral range, the second due to Compton scatter¬ 
ing peaking in the MeV-GeV spectral range. The y-weak AGN, on the other hand, are 
the less luminous Seyfert and radio galaxies (being spiral and ellipticals, respectively), 
where the spectra extends from radio frequencies to the soft y-rays (100s of keV). These 
spectra have approximately equal power per logarithmic frequency interval, but show 
many spectral features, in particular, three broad band spectral peaks at infrared, UV, 
and hard X-ray frequencies. The fast time variability of the UV to soft y-ray spectra 
indicates that this emission is ’’nuclear” radiation originating in the accretion flow very 
close to the black hole. The slowly varying broad band infrared spectral peak originates 
at larger distances, most likely being due to nuclear radiation being reprocessed by dust. 
The nuclear spectra has been studied for more than 20 years, yet little is known about 
the physical conditions close to the super massive black holes in AGN. 

Both Seyfert (normally being radio-weak) and radio galaxies comes in two types based 
on optical classification. Seyfert 1 and the broad line radio galaxies show Doppler- 
broadened broad (up to 10 000 km/s) and narrow (up to perhaps 1000 km/s) emission 
lines, while the Seyfert 2 galaxies and the narrow line radio galaxies show only the narrow 
emission lines. The broad lines originate in fast-moving photo-ionized 10^ K gas clouds 
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within 1 parsec of the central source, while the narrow lines originate in slower moving 
clouds at kiloparsec distances. Aspects of radio galaxies are discussed in the chapter by 
Malcolm Longair in this volume. Here, we will focus on the Seyfert galaxies only. 

In the unified scheme for Seyfert galaxies (e.g. Antonucci 1993), we are directly viewing 
the central X-ray source and broad line region in Seyfert 1 galaxies, but our line-of-sight 
passes through an optically opaque molecular torus in Seyfert 2 galaxies obscuring both 
the central source and the broad line region but not the narrow line region. The photo- 
ionized narrow line gas lies in an ionization cone whose half-opening angle is determined 
by the geometry of the obscuring torus. In some Seyfert 2 galaxies we still see weak 
polarized broad lines being reflected by electrons in the ionization cone. In standard 
accretion disk scenarios for the central X-ray source, the unified models imply that our 
viewing angle (i.e. the angle between the disk normal and the line-of-sight) is less than the 
half-opening angle of the ionization cone for the case of Seyfert Is. Half-opening angles 
inferred from observations are typically 30-40 degrees, which means that in Seyfert Is we 
are viewing the accretion disk from directions that are closer to face-on than to edge-on. 


3. X-ray and 7 -ray spectra of Seyfert galaxies 

The satellites (launch year and energy range in parenthesis) ROSAT (June 1990; 0.1-2 
keV), ASCA (Feb 1993; 0.1-10 keV), Ginga (Feb 1987 - Nov 1991; 2-40 keV) , SAX 
(April 1996; 0.2-200 keV), XTE (Dec 1995; 2-200 keV), and CGRO (April 1991; 50 keV- 
30 GeV) have allowed broad band X/ 7 -ray studies of some 30 Seyfert galaxies in the 
1990s. However, as AGN varies in the X-rays on time-scales down to days or even hours 
simultaneous measurements are necessary. Such simultaneous observations unfortunately 
exist for only a handful objects, but the number is increasing. The launch of INTEGRAL 
(2 keV-10 MeV) in the year 2001 will finally provide a large number of truly simultaneous 
X/ 7 -spectra of AGN. 

By combining observations in different energy bands, a picture has emerged where the 
overall shapes of the X-ray and 7 -ray spectra of different Seyfert 1 galaxies are similar. 
The only Seyfert 1 galaxy for which the broad band X- and 7 -ray spectrum is well de¬ 
termined is IG 4329A (Madejski et al. 1995) being the second brightest Seyfert galaxy 
in the 2-10 keV range. Simultaneous ROSAT and CGRO OSSE, as well as nonsimulta- 
neous Ginga data are shown in Figure [left panel). The observations were modelled by 
Madejski et al. (1995), Zdziarski et al. (1994), and Magdziarz & Zdziarski (1995). The 
spectrum consists of two components: 1) an intrinsic power law of slope Oint ~ 0.9 with 
an exponential cutoff energy of about 300 keV (dashed curve in Fig. |^), and 2) a reflection 
component caused by cold reflecting matter subtending a solid angle ~ 1 — 27r as viewed 
from the X-ray source (dotted curve in Fig. |^. The reflection component consists of a 
fluorescent Fe line at 6.4 keV and a broad peak at about 30 keV causing the overall 2-18 
keV spectral slope to be considerably flatter, a 2 -i 8 ~ 0.7 and the apparent cutoff energy 
to be smaller, « 50 keV. A small column depth of neutral material about Nn = 3 x 10^^ 
cm“^ along the line of sight causes absorption below 2 keV. 

The average X/ 7 -spectrum for four Seyfert Is using nonsimultaneous Ginga and OSSE 
data shows similar properties as the spectrum for IG 4329A (Zdziarski et al. 1995). The 
same holds for the nonsimultaneous average EXOSAT/OSSE spectrum for seven Seyfert 
Is (Gondek et al. 1996). 

Larger samples exist within narrower spectral ranges using a single experiment. These 
samples show spectra consistent with the smaller samples above. For example, 60 Ginga 
spectra of 27 Seyferts show aint, 2 -i 8 ~ 0.95 and the overall 02-18 ~ 0.73, each having a 
dispersion of about 0.15 (Nandra & Pounds 1994). 


R. Svensson: X-rays and gamma rays from AGN 


5 



Figure 1. Left panel shows the observed broad band spectrum (crosses) of IC 4329A from 
ROSAT, Ginga, and CGRO OSSE (Madejski et al. 1995). The dashed curve shows the intrinsic 
model spectrum incident on both the cold matter and the observer, the dotted curve is the 
reflected component, and the solid curve is the observed sum. From Magdziarz & Zdziarski 
(1995). Right panel shows the observed broad band spectrum of NGC 4151 from Ginga, and 
GGRO OSSE (crosses) and from EXOSAT (circled crosses). The EXOSAT data are from April 
1984, the highest and lowest Ginga data from May 1990 and May 1987, respectively, and the 
highest and lowest GGRO OSSE data from May and April 1993, respectively. Solid, dashed, 
and dotted curves show hts to the data using a simple thermal Comptonization model. From 
Zdziarski, Johnson & Magdziarz (1996). 


The brightest Seyfert galaxy, the Seyfert 1.5 NGC 4151, once considered a prototype 
for Seyfert Is, has now turned out to have unusua l X-ray properties. Figure (right 
panel) from Zdziarski, Johnson fc Magdziarz (1996) shows the extreme spectral states of 
NGC 4151. In the 2-18 keV Ginga-iange, the spectral index, a 2 -i 8 , varies in the range 
0.3-0.8 instead of having the canonical Seyfert value of 0.7. The harder the spectrum 
is, the weaker is the 2-18 keV flux, and the spectrum seems to pivot around 100 keV. 
And indeed, OSSE observations during 1991-1994 show that the OSSE flux is essentially 
constant. NGC 4151 shows no or a smaller reflection component than normal Seyfert Is. 
This causes the deduced intrinsic power law component to have a smaller cut off energy 
than the standard Seyfert 1, although the total OSSE spectra are indistinguishable. NGC 
4151 also have a larger column depth A^h ~ 10^^ cm“^ of absorbing gas along the line of 
sight compared to typical Seyfert Is. This causes absorption below about 6 keV. Optical 
obervations of an ionization cone indicates that we are not viewing NGC 4151 face-on 
but rather at a viewing angle of 65 degrees. NGC 4151 is clearly a freak object which 
requires special consideration to fit into the unified scheme of Seyfert galaxies. The 
standard interpretation is that we have a direct view of the central X-ray source, but it 
has also been proposed (Poutanen et al. 1996) that the central X-ray source is obscured 
and that the observed X-ray spectrum is a standard Seyfert 1 spectrum scattered by 
about 65 degrees into the line of sight. Compton recoil then decreases the spectral cut 
off energy. 

The observed obscuring column depths in Seyfert 2 galaxies lies in the range IVh ~ 
2 X 10^^ - 5 X 10^^ cm“^ (Smith & Done 1996). An extreme case is NGC 4945 (Done, 
Madejski & Smith 1996) whose broad band X-ray spectrum is shown in Figure |^. Here, 
the absorbing column is very large, A^h ~ 4 x lO^'^ cm“^, causing absorption of the direct 
spectrum below 20 keV. Below 10 keV the spectrum is dominated by a weak component 
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Figure 2. The broad band spectrum (crosses) of NGC 4945 from ASCA, Ginga, and CGRO 
OSSE. The dotted eurve is the direct power law spectrum being absorbed below 20 keV by a 
column depth A^h ~ 4 x 10^^ cm“^. The dash-dotted eurve shows the component scattered into 
the line-of-sight by electrons in the ionization cone. The effects of a reflection component is 
included. From Done, Madejski & Smith (1996). 

that has been scattered into the line-of-sight by electrons in the ionization cone. Although 
weak in the 2-10 keV range, NGC 4945 is the second brightest known Seyfert galaxy in 
the sky above 50 keV. In the Seyfert 2 galaxy, NGC 1068, the absorbing column is even 
larger and the scattered component dominates at least up to 20 keV (Koyama et al. 
1989). 

4. Outstanding questions regarding the X-ray emitting regions in 
Seyfert 1 galaxies 

The obvious question is of course which radiation process generates the X-ray con¬ 
tinuum. It is commonly believed that Compton scattering by energetic electrons, either 
mildly relativistic thermal electrons or highly relativistic nonthermal electrons, Compton 
scatter soft UV-photons into the X-ray range. 

Then there is the question of geometry, i.e. the spatial distribution of hot X-ray 
generating electrons and of cold (T < 10® K) reflecting matter. Possible geometries 
include a cold slab surrounded either by plane parallel coronal slabs of hot electrons, or 
by coronal patches (active regions) of unknown geometries. 

Finally, there is the question whether it is thermal or nonthermal electrons or both 
that account for the X-ray emission. This question might be answered if the spectrum 
above a few hundred keV was known with certainty, but the best available spectrum 
of the brightest typical Seyfert 1, IC 4329A, does not have sufficient signal to noise at 
such energies. The spectrum of IC 4329A can be fit with both thermal and nonthermal 
models (Zdziarski et al. 1994). 

The nonthermal models studied in the 1980s have several attractive features (for re¬ 
view, see Svensson 1994). In particular, it was predicted already in 1985 (Svensson 1986) 
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that the pair cascades in the nonthermal models give rise to an X-ray spectrum with 
aint, 2 -i 8 ~ 0.9— 1.0 in contradiction to the overall slope of a 2 -i 8 ~ 0.7 but in agreement 
with the intrinsic spectral spectral slope later found by Ginga (e.g. Nandra & Pounds 
1994). 

Due to the non-detections by COMPTEL and the high energy cutoffs indicated by 
OSSE, attention is now focused on the thermal models. 


5. Geometry of the X-ray source 

In principle, the central X-ray source could have spherical symmetry with the source of 
UV-photons being small cold clouds uniformly distributed throughout the X-ray source. 
In this case, radiation models invoking spherical symmetry would be applicable (e.g., 
Ghisellini & Haardt 1994, Skibo et al. 1995; Pietrini & Krolik 1995). In some models, 
the central X-ray source is quasi-spherical but the UV-source is anisotropic. Here, the X- 
ray source can either be a hot two-temperature SLE-accretion disk or a two-temperature 
ion torus. The UV-source may then be the standard cold Shakura-Sunyaev disk located 
beyond some transition radius. Another proposed alternative is that the disk is hot at 
larger radii but that in the innermost region the cooling times become sufficiently short 
for the disk to collapse into a cold UV-emitting disk state. The resulting anisotropy 
of the UV-photons breaks the spherical symmetry in both of these scenarios. Then 
detailed radiative transfer calculations of how the UV-photons penetrate the hot disk are 
required in order to determine the emerging X-ray spectrum of Comptonized UV-photons 
at different viewing angles. This problem has not yet been solved. 

The most commonly used scenario (mainly because of its simplicity) is the two-phase 
disk-corona model (e.g., Haardt & Maraschi 1991, 1993) in which a hot X-ray emitting 
corona is located above the cold UV-emitting disk of the canonical black hole model for 
AGNs. The power law X-ray spectrum is generated by thermal Comptonization of the 
soft UV-radiation. About half of the X-rays enters and is reprocessed by the cold disk, 
emerging mostly as black body disk radiation in the UV. Haardt & Maraschi (1991) 
emphasized the coupling between the two phases due to the reprocessing, as the soft disk 
photons influence the cooling of the corona. They showed that nearly all power must 
be dissipated in the corona in order to have aint ~ 0.9 — 1. A consequence of this is 
that the soft disk luminosity, Lg, is of the same order as the hard X-ray luminosity, Lh- 
The disk-corona scenario is highly anisotropic as the UV-photons enter the corona from 
below only. 

Observations show that Lg may be several times larger than Lh, in contradiction to 
the prediction of the uniform two-phase disk-corona model. This led Haardt, Maraschi, 
& Ghisellini (1994) to propose a patchy disk-corona model, where the corona consists of 
several localized active regions on the surface of the disk. Internal disk dissipation results 
in UV-radiation that leaves the disk without entering the active regions, thus enhancing 
the observed Lg relative to Lh- The patchy corona model is also highly anisotropic as 
the UV-radiation enters the active regions from below and through the sides. The active 
regions could also be elevated above the disk. 

The full radiative transfer and Comptonization problem in the these geometries has 
recently been solved (Haardt & Maraschi 1993, Stern et al. 1995b, Poutanen & Svensson 
1996). We now turn to discuss the methods and the results. 
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6. Radiative transfer/Comptonization 

Two different methods have been used to solve the full radiative transfer/Comptoniza- 
tion problem in mildly relativistic thermal plasmas accounting for energy and pair balance 
as well as reprocessing by the cold disk (including angular anisotropy and Klein-Nishina 
effects). 

The first method is based on the Non-Linear Monte Carlo (NLMC) method developed 
by Stern (1985, 1988) and described in detail in Stern et al. (1995a). The Monte Carlo 
particles (particles and photons) make up the background so the Monte Carlo particles 
can interact with themselves (e.g. photon-photon pair production) causing nonlinearity. 
Any geometry can be treated including 2D and 3D geometries, but calculations so far 
(Stern et al 1995b) have been limited to coronal slabs (ID) or active regions (2D) in 
the shape of hemispheres or spheres at different elevations above the cold disc. Another 
advantage is the possibility to divide the region into several zones in order to study the 
inhomogeneous distributions of, e.g., temperature and pair density. The drawback is that 
each run takes a few hours on a Sun IPX. 

The second method is based on the iterative scattering method (ISM), where the ra¬ 
diative transfer is exactly solved for each scattering order separately (e.g. Sunyaev & 
Titarchuk 1985, Haardt 1994). The code and its testing are briefly described in Poutanen 
and Svensson (1996, see also Poutanen 1994; Poutanen & Vilhu 1993). It was shown that 
the results of the NLMC and the ISM code are in excellent agreement. The ISM code is a 
ID code but it can also treat quasi-2D radiative transfer in cylinders/pill boxes and hemi¬ 
spheres. The full Compton scattering matrix is used allowing to solve polarized radiative 
transfer in thermal relativistic plasmas. Fully angular dependent, polarized Compton re¬ 
flection is implemented using a Green’s matrix (Poutanen, Nagendra & Svensson 1996). 
The advantage of the ISM code is that it takes 10 minutes or less on a Sun IPX. 


7 . The Comptonization solution 

The two codes have been used to study radiative transfer/Comptonization in pure 
pair coronae in energy and pair balance. For coronae of a given geometry and in energy 
balance, there exists a unique Te — tt relation, where Te is the volume-averaged coronal 
temperature and tt is a characteristic Thomson scattering optical depth of the coronal 
region. For slabs, ry is the vertical Thomson scattering optical depth. For hemispheres 
and for spheres, the radial tt is averaged over 27r solid angle for hemispheres, and Att 
solid angle for spheres. 

In Figure [^a, this relation is shown for different geometries. The results for slabs 
are shown by rectangles, for hemispheres located on the disk surface by hemispheres, 
for surface spheres located on the disk surface by underlined spheres, and for spheres 
located at heights, 0.5h, Ih, and 2h, by spheres, diamonds, and triangles, respectively 
[h is the radius of the sphere). The results for active regions are connected by dotted 
curves. For comparison we also show the slab results from Stern et al. 1995b using the 
ISM code {dashed curve). Each curve is characterized by an almost constant generalized 
Kompaneets parameter, y = rT(l -I- rT)(40 -|- 160^), where 0 = kTe/meC^, and where 
y ~ 0.49 for slabs, ~ 2.0 for hemispheres, and ~ 2.3, 3.7, 5.9, and 7.8 for the spheres. The 
larger the soft photon starvation (i.e. the fewer the number of reprocessed soft photons 
reentering the coronal region), the larger is y. There is very good agreement between the 
slab results from the NLMC and the ISM codes, which tests the accuracy of both codes 
and methods. 


R. Svensson: X-rays and gamma rays from AGN 


9 



Figure 3. Dimensionless volume-averaged temperature, 0 = kTejm^c?, vs. Thomson scatter¬ 
ing optical depth, tt, in panel (a), and vs. dissipation compactness, £diss = {Ldiss/h) (aT/meC^) 
in panel (b), for a steady X-ray emitting plasma region in pair and energy balance on or above 
a cold disk surface of temperature, kT^h = 5 eV. The plasma Compton scatters reprocessed soft 
black body photons from the cold disk surface. Solid rectangles and dashed curve show results 
from NLMC code and ISM code, respectively, for the case of a plane-parallel slab corona. Re¬ 
sults using the NLMC code for individual active pair regions are shown for hemispheres located 
on the disk surface; surface spheres also located on the surface {underlined spheres); spheres 
located at a height of 0.5/i {spheres), Ih {diamond), and 2h {triangle), where h is the radius 
of the sphere. The results for each type of active region are connected by dotted curves. The 
dash-dotted and dash-dot-dot-dotted curves in panel (b) show the critical compactness as func¬ 
tion of 0 above which thermalization by Moller and Bhabha scattering is not achieved for the 
cases of pair slabs and surface spheres, respectively. From Stern et al. (1997). 


8. Pair balance and the compactness 

Solving the pair balance for the combinations of (0, tt) obtained in § ^ gives a unique 
dissipation compactness, ^diss (see Ghisellini & Haardt 1994 for a discussion). Here, the 
local dissipation compactness, £diss = {Ldiss/h){a'T/meC^), characterizes the dissipation 
with Ldiss being the power providing uniform heating in a cubic volume of size h in the 
case of a slab of height h, or in the whole volume in the case of an active region of size h. 
The Figure IHd shows the 0 vs. £diss relation obtained with the NLMC code for different 
geometries and for fcTbb = 5 eV. The slab results with the ISM codes are also shown 
by the dashed curve. The parameter space to the right of respective curve is forbidden 
as pair balance cannot be achieved, and the parameter space to the left would contain 
solutions where the background coronal plasma dominates over the pairs, i.e. “pair free” 
solutions (e.g. Svensson 1984, HM93). Note that the pair temperatures are locked in the 
range 0 ^ 0.15 — 1 (80 - 500 keV) for the considered compactnesses. The pairs thus act 
as a “thermostat” (e.g., Svensson 1984) 

From Figure l^b we find that at a given 0 the active regions have a larger £diss than the 
slabs. This is due to the longer escape times in slabs. In a local cubic volume in slabs, 
the four vertical surfaces act effectively as “reflecting” surfaces increasing the radiation 
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energy density in the volume as compared to the active regions where the radiation 
escapes freely through the surface of the volume. 


9. Thermalization 

The question arises whether the electrons can thermalize or not for the conditions, 
0, tt, and ^diss, in Figure 0. Energy exchange and thermalization through Moller 


(e 




and Bhabha (e^e ) scattering compete with various loss mechanisms, with 


Compton losses being the most important for our conditions. The thermalization is slow¬ 
est and the Compton losses largest for the higher energy particles in the Maxwellian tail. 
Ghisellini, Haardt fc Fabian (1993)| compared average Maxwellian time scales to find the 
conditions when thermalization is not achieved. Instead we use the detailed simulations 


by Dermer & Liang (1989) (their Fig. 8) to find the critical compactness above which the 
deviation at the Maxwellian mean energy is more than a factor e « 2.7. The dash-dotted 
and dash-dot-dot-dotted curves in Figure ^ show this critical compactness for slabs and 


for surface spheres, respectively. In agreement with Ghisellini, Haardt & Fabian (1993) 


we find that Moller and Bhabha scattering cannot compete with Compton losses in our 


pair slabs and active regions. On the other hand, Ghisellini, Guilbert & Svensson (1988) 


and Ghisellini fc Svensson (1990) found that cyclo/synchrotron self-absorption acts as a 
very efficient thermalizing mechanism as long as the magnetic energy density dominates 
the radiation energy density, a condition expected in the corona and in the magnetic 
flares on the surface of accretion disks. 


10. Anisotropy effects 

Figure ^ shows the “face-on” emerging spectrum for a hemisphere with £diss = 20 and 
fcTbb = 5 eV with the reflected component subtracted. The solid curve is the sum of 
the reprocessed soft black body spectrum and the Comptonized spectrum. The “face- 
on” spectrum is averaged over viewing angles 0.6 < cos9 < 1. The numbered dotted 
curves show spectral profiles of the different scattering orders. As noted by Haardt 
(1993) and HM93, the first scattering order is significantly reduced in face-on directions 
due to anisotropic Compton scattering when 0 is mildly relativistic. This deficiency at 
low energies causes an anisotropy break close to the peak of the second scattering order. 
Below the anisotropy break, the spectrum is a power law with an anisotropy slope, Oaniso, 
that is harder than the standard isotropy slope, Oiso, above the break. It is the standard 
isotropy slope, Oiso) that has been fitted with several analytical expressions that are 
functions of 0 and tt (e.g., Zdziarski 1985, Titarchuk 1994). “Edge-on” spectra have a 
much weaker reflection component and no anisotropy break. 


11. Spectra from active pair regions: hemispheres 

Figure || shows emerging “face-one” and “edge-on” spectra from hemispheres at dif¬ 
ferent £diss for fcTbb = 5 eV. We have Oaniso ~ 0.77 while aiso ~ 1.09. The anisotropy 
break is therefore about 0.3. The anisotropy break moves through the 2-18 keV range 
for .^diss ~ 10-100, and thus 02-18 is smaller than unity for ^diss less than 100. At this 
low fcTbb the first Compton order hardly extends into the 2-18 keV range and there is 
little discreteness effect on the behavior of 02 - 18 - 
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log e 

Figure 4. Emerging spectrum, eFl, where is the energy flux (arbitrary units) and 
e = hvfmeC? from hemisphere with compactness, faiss = 20, and cold disk temperature, kT^h 
= 5 eV. The spectrum is averaged over viewing angles 0.6 < cos 9 < 1. The solid curve shows 
the sum of the Comptonized spectrum from the hemisphere and the reprocessed black body 
spectrum (the reflection component is not included). The numbered dotted curves show the 
different scattering orders. Note that the first scattering order is significantly reduced due to 
anisotropic Compton scattering. This deficiency causes an anisotropy break close to the peak of 
the second scattering order. Below the anisotropy break, the spectrum is a power law with an 
anisotropy slope that is harder than the standard isotropy slope above the break. From Stern 
et al. (1997). 

12. Diagnostics Using Compactness and Ginga Slopes 

The least square overall spectral slope, a 2 _i 8 , for the 2-18 keV range were determined 
and are displayed in Figure || as a function of the dissipation compactness, f^diss: for 
different geometries. The right panel of Figure ^ shows the observed distribution of 
02-18 for Ginga spectra from 27 Seyfert galaxies (Nandra & Pounds 1994). The crosses 
represent 17 Seyfert galaxies that have both been observed by Ginga and have known 
estimates of their X-ray time variability (and thus lower limits of their compactnesses). 
The true crosses may lie to the right of the plotted ones. 

One sees that the observations are more consistent with active surface regions (such as 
hemispheres or surface spheres) than with slabs. Active regions produce spectra covering 
the observed ranges of 02-18 (~ 0-4 — 0.9) and cutoff energies (~ 2kTg) for the observed 
range of compactnesses. 


13. Conclusions 

Have we been able to throw some light upon the outstanding questions discussed in 

Regarding the radiative transfer/Comptonization, the progress has been rapid and it 
is now possible to use full radiative transfer simulations when fitting X-ray observations, 
e.g. with XSPEC, at least for the simplest geometries (e.g., Poutanen et al. 1996 for the 
case of NGC 4151). As anisotropic effects are very important, i.e. as spectral shapes 
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Figure 5. Emerging “face-on” and “edge-on” spectra, tF^, where FI is the energy flux (ar¬ 
bitrary units) and e = hvjmeC?, from hemispheres of different compactnesses, idiss and for 
cold disk temperature, kThh = 5 eV. The “face-on” spectra are averaged over viewing angles 
0.6 < cos 6 < 1, and the “edge-on” spectra over viewing angles 0 < cosO < 0.5. The solid curves 
show the total spectrum which is the sum of the Comptonized spectrum from the hemisphere 
itself {dotted eurves), the reprocessed black body spectrum, and the reflection component. Ver¬ 
tical dashed lines show the 2-18 keV spectral range. The anisotropy break in the face-on spectra 
moves to lower photon energies as idiaa increases. From Stern et al. (1995b, 1997). 


depend strongly on viewing angle, it will be possible to set constraints on the viewing 
angle. The ISM code also gives detailed predictions for the spectral X/y-ray polarization, 
which may provide crucial constraints once observations become available. 

The spectra also depend on the geometry of the coronal regions, so observed spectra 
can be used as diagnostics of the geometry. Presently, it seems that active regions are 
favoured over homogeneous slab coronae. 

There is still no answer to the question whether the radiating electrons (and positrons) 
are thermal or nonthermal. High quality spectra above 200 keV are needed, and the 
question may not be settled before the launch of INTEGRAL in the year 2001. 
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sions. Boris Stern’s research visits at the Stockholm Observatory were financed by the 
Swedish Institute, the Swedish Royal Academy of Sciences, and by a NORDITA Nordic 
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